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ABSTRACT 

Statistics of lensed arcs in clusters of galaxies serve as a powerful probe of both the non-sphericity and the 
inner slope of dark matter halos. We develop a semi-analytic method to compute the number of arcs in triaxial 
dark matter halos. This combines the lensing cross section from the Monte Carlo ray-tracing simulations, and 
the probability distribution function (PDF) of the axis ratios evaluated from cosmological A^-body simulations. 
This approach enables one to incorporate both asymmetries in the projected mass density and elongations along 
the line-of-sight analytically, for the first time in cosmological lensed arc statistics. As expected, triaxial dark 
matter halos significantly increase the number of arcs relative to spherical models; the difference amounts to 
more than one order of magnitude while the value of enhancement depends on the specific properties of density 
profiles. Then we compare our theoretical predictions with the observed number of arcs from 38 X-ray selected 
clusters. In contrast to the previous claims, our triaxial dark matter halos with inner density profile p oc r" 1 5 in 
a Lambda-dominated cold dark matter (CDM) universe reproduces well the observation. Since both the central 
mass concentration and triaxial axis ratios (minor to major axis ratio ~ 0.5) required to account for the observed 
data are consistent with cosmological N-body simulations, our result may be interpreted to lend strong support for 
the collisionless CDM paradigm at the mass scale of clusters. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general — gravitational lensing 

1. INTRODUCTION 

The discovery of a lensed arc in a rich cluster A370 (Lynds & Petrosian 1986; Soucail et al. 1987) opened a direct window to 
probe the dark mass distribution in clusters of galaxies. Since gravitational lensing phenomena are solely dictated by intervening 
mass distributions, they are not biased by the luminous objects unlike other conventional observations. Indeed, previous work (Wu 
& Hammer 1993; Miralda-Escude 1993a, 1995, 2002; Bartelmann 1996; Hattori, Watanabe, & Yamashita 1997b; Molikawa et al. 
1999; Williams, Navarro, & Bartelmann 1999; Meneghetti et al. 2001; Molikawa & Hattori 2001; Oguri, Taruya, & Suto 2001; 
Sand, Treu, & Ellis 2002; Gavazzi et al. 2003) showed that the number, shapes, and positions of lensed arcs are sensitive to the 
mass distribution of clusters. For instance, Oguri et al. (2001) calculated the number of arcs using the generalization of the universal 
density profile proposed by Navarro, Frenk, & White (1996, 1997) and pointed out that it is extremely sensitive to the inner slope and 
the concentration parameter of the density profile; the number of arcs changes by more than an order of magnitude among different 
models that are of cosmological interest. Therefore lensing arc surveys provide an important probe of density profiles of clusters in 
a complementary manner to the statistics of wide-separation lensed quasars (Maoz et al. 1997; Keeton & Madau 2001; Oguri 2002b, 
2003; Li & Ostriker 2003). 

While most previous studies of lensed arcs have aimed at constraining the cosmological parameters (Wu & Mao 1996; Bartelmann 
et al. 1998; Cooray 1999; Sereno 2002; Golse, Kneib, & Soucail 2002; Bartelmann et al. 2003), we rather focus on extracting 
information of the density profiles of dark matter halos. Thus we assume a Lambda-dominated cold dark matter (CDM) model that is 
consistent with the recent Wilkinson Microwave Anisotropy Probe (WMAP) result (Spergel et al. 2003); the matter density parameter 
flo = 0.3, the dimensionless cosmological constant Ao = 0.7, the mass fluctuation amplitude a% = 0.9, the Hubble constant in units of 
100kms~ 1 Mpc" 1 , h = 0.7. In fact, arc statistics depend on the assumed set of cosmological parameters in two ways; directly through 
the geometry of the universe and somewhat indirectly through properties of density profiles which also depend on the cosmology. 
For instance, Bartelmann et al. (1998) found that the numbers of arcs significantly change among different cosmological models, and 
concluded that only open CDM models can reproduce the high frequency of observed arcs. Oguri et al. (2001) showed, however, 
that the result largely comes from the larger concentration parameter of halo profiles in the open CDM model than in the Lambda- 
dominated CDM model. Thus this may be more related to the small-scale behavior of the CDM model than the "global" effect of the 
cosmological constant. 

The values of cosmological parameters are determined fairly accurately now, thus our primary interest here is to confront the 
density profiles of dark matter halos with the arc statistics, and thereby we would like even to test the collisionless CDM paradigm. 
For this purpose, a non-spherical description for the lensing halos is the most essential since cross sections for arcs are quite sensitive 
to the non-sphericity of mass distribution (e.g., Bartelmann, Steinmetz, & Weiss 1995; Bartelmann 1995; Meneghetti et al. 2001; 
Oguri 2002a). Indeed, previous analytic models adopting spherical lens models failed to reproduce the observed high frequency of 
arcs (Hattori et al. 1997b; Molikawa et al. 1999). Because of the lack of a realistic analytical model for non-spherical lens, however, 
one had to resort to N-body simulations to take account of non-spherical effects on the arcs statistics (e.g., Bartelmann & Weiss 
1994). Nevertheless it is quite demanding for those numerical simulations to resolve the central part of the gravitational potential of 
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lensing halos while keeping the reasonable number of those objects sufficient for statistical discussion. This is why a complementary 
(semi-) analytical approach to the arc statistics is highly desired. 

Recently, new methods to constrain the mass profile of individual clusters also have been developed (Smith et al. 2001; Sand et al. 
2002; Clowe & Schneider 2001, 2002; Gavazzi et al. 2003). Although such methods can measure mass distributions of individual 
clusters precisely, it may suffer from the special selection function and the scatter around the mean mass distribution. For instances, 
analysis of clusters only with giant arcs may result in more elongated clusters than average because JS02 showed that triaxial axis 
ratios have fairly broad distributions. Therefore it is of great importance to study statistics of lensed arcs which allow us to obtain 
information on the mean profile. 

In this paper, we develop and study in detail, for the first time, such an analytical model of the non-spherical lensing objects for 
the arc statistics. Specifically we adopt the triaxial description of dark matter halos proposed by Jing & Suto (2002, hereafter JS02). 
They have presented detailed triaxial modeling of halo density profiles, which enables us to incorporate the asymmetry of dark matter 
halos statistically and systematically. We first compute the lensing cross sections for arcs on the basis of the Monte Carlo simulations 
following Oguri (2002a). Then we make systematic predictions of the number of arcs by averaging the cross sections over the 
probability distribution functions (PDFs) of the axis ratios and the concentration parameters and assuming the random orientation of 
the dark halos along the line-of-sight of the observer. Those theoretical predictions are compared with the number of observed arcs in 
a sample of 38 X-ray selected clusters compiled by Luppino et al. (1999). We pay particular attention to several selection functions 
of clusters and arcs which may systematically affect our results (e.g., Wambsganss, Bode, & Ostriker 2003). 

The plan of this paper is as follows. In §2, we briefly summarize the triaxial modeling proposed by JS02. We present several key 
results of gravitational lensing in triaxial dark matter halos in §3. The method to predict the number of arcs is described in detail in 
§4, and comparison with observations is discussed in §5. Finally, we discuss several implications of our results in §6, and summarize 
the conclusion in §7. 



2. DESCRIPTION OF TRIAXIAL DARK MATTER HALOS 

In this section, we briefly summarize the triaxial model of dark matter halos proposed by JS02. They obtained the detailed triaxial 
modeling on the basis of their high-resolution individual halo simulations as well as large-scale cosmological simulations. Most 
importantly, they provided a series of useful fitting formulae for mass- and redshift-dependence and the PDFs of the axis ratio and 
the concentration parameter. Such detailed and quantitative modeling enables us to incorporate the non-sphericity of dark matter 
halos in a reliable manner. 



2.1. Coordinate Systems 

We introduce two Cartesian coordinate systems, x = (x, y, z) and x' = (x 1 , y' , z')> which represent respectively the principal coordinate 
system of the triaxial dark halo and the observer's coordinate system. The origins of both coordinate systems are set at the center of 
the halo. It is assumed that the z'-axis runs along the line-of-sight direction of the observer, and that the z-axis lies along the major 
principal axis. In general, the relative orientation between the two coordinate systems can be specified by the three Euler angles. 
However, in our case, it is only the line-of-sight direction that is fixed while the rotation angle of the x'-y' plane relative to x-y plane 
is arbitrary, and thus we may need only two angles to specify the relative orientation of the two coordinate systems. Here we make a 
choice of x'-axis lying in the x-y plane. Then the relative orientation of the two coordinate systems can be expressed in terms of the 
line-of-sight direction in the halo principal coordinate system. 

Let (9, <j>) be the polar coordinates of the line-of-sight direction in the incoordinate system. Then the relation between the two 
coordinate systems can be expressed in terms of the rotation matrix A (Binney 1985) as 

x=Ax', (1) 

where 

-cos</>cos# coscf>sm6 \ 

-sin(/>cos6> sin^sin^ . (2) 
sin8 cos 8 ) 

Figure 1 represents the relative orientation between the observer's coordinate system and the halo principal coordinate system. 

2.2. Density Profile of Triaxial Dark Matter Halos 
We adopt the following density profiles of triaxial dark matter halos proposed by JS02: 

<Wcrit(z) 

(R/R r(l + R/Ro) 

where 




(W = ,„/„ . n/n ,3- a ' ( 3 ) 



*W(^44) (a< b < c) . (4) 

The precise value of the inner slope, a, is still controversial, but almost all the N-body simulations based on the collisionless CDM 
scenario indicate values between 1 and 1.5 (Navarro et al. 1996, 1997; Moore et al. 1999; Jing & Suto 2000; Fukushige & Makino 
2001, 2003; Power et al. 2003). Thus we consider both a = 1 and a = 1 .5 below so as to cover a possible range of the CDM predictions 
proposed so far. 
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FIG. 1. — The orientations of the coordinate systems. The Cartesian axes (x,y,z) represent the halo principal coordinate system while the axes (x' ,y' ,z') stand for 
the observers coordinate system with z'-axis aligned with the line-of-sight direction. The .v'-axis lies in the x-y plane. The angle (£>, <j>) represent the polar angle of 
the line-of-sight direction in the (x,y, z)-coordinate system. 



JS02 defined the concentration parameter in the triaxial model as 



where R e is chosen so that the mean density within the ellipsoid of the major axis radius R e is A e f2(z)p cr i t (z) with 2 

0.75 



A, = 5A V 



ab J 



(5) 



(6) 



Here il(z) and p a i t (z) denote the matter density parameter and the critical density of universe at redshift z, respectively, and A v ; r (z) 
denotes the overdensity of objects virialized at z whose approximate expression is found, e.g., in Oguri et al. (2001). 
Then the characteristic density <5 ce in equation (3) is written in terms of the concentration parameter c e as 

r A e ft(z) c] 



m(c e ) ' 



where m(c e ) is 



„3-Q 



m(c e ) ■ 



■ %F\ (3 — a, 3 — a; 4 — a; — c e ) , 



3-a' 

with 2F1 (a,b;c;x) being the hypergeometric function. For a = 1 and 1.5, equation (8) simply reduces to 

c. 



ln(l+c e )- 



m(c e ) = 



l+c e 



(a=l), 



2\n{JT e + JY+c e )-2j (a =1.5) 
V l+c e 



(7) 



(8) 



(9) 



Since R e is empirically related to the (spherical) virial radius r wu as R e /r vn ~ 0.45 (JS02), the scaling radius in the triaxial model, Rq, 
for a halo of a mass M W1I is given as 



^, = 0.45—: 



0.45 



3M V 



47rA vir f2(z)p cl - it (z) 



1/3 



(10) 



Since we do not know the properties of the density profile of an individual lensing halo, our prediction for the number of arcs is 
necessarily statistical in a sense that it should be made after averaging over appropriate PDFs of the properties of halos. For this 
purpose, we adopt the PDFs that JS02 empirically derived from their simulations. For the axis ratios, they are given as 

{( fl / C )(M vn ./M*) 007 [^]° 7 - 0.54} 2 " 



p(a/c) ■■ 



1 

2ttx0.113 



exp 



2(0.113) 2 



M vh 
AT, 



0.07[f2(z)]° 



(ID 



2 Note that our definitions of A v j r and A e are slightly different from those of JS02; A v j r (JS02) = n(z)A v j,, and A e (JS02) = f2(z)A E . Of course this does not change 
the definition of R e . 
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and 



p(a/b\a/c)-- 



2(l-max(a/c,0.5)) 



2a/ b — 1 -max(a/c,0.5) 
1 -max(a/c,0.5) 



(12) 



for a/b > max(a/c,0.5), and p(a/b\a/c) = otherwise (JS02). Here M* is the characteristic nonlinear mass so that the rms top-hat 
smoothed overdensity at that mass scale is 1 .68. For the concentration parameter, we adopt 

1 I" (lnc e -lnc e ) 21 

P( c e) = TT ex P 



2(0.3) 2 



/27T x 0.3 

where the fit to the median concentration parameter c e for a = 1 is given as 3 : 

! 11.1 I 

c e = 1.35exp 



1 



(a/c)(M vil /M*)™W(z)f 7 




(13) 



(14) 



with z c being the collapse redshift of the halo of mass M v ; r (JS02). In the case of a = 1, we simply use the above expression, and 
for a = 1.5, we use the relation c e (a = 1.5) = 0.5c e (a = 1) (Keeton & Madau 2001; JS02). JS02 estimated A e = 1.1 in the Lambda- 
dominated CDM model, but this value is likely to be dependent on the underlying cosmology to some extent. As we stressed in 
Introduction, however, we do not intend to survey the cosmological parameters but rather focus on the effects of the properties of the 
lensing halos. Therefore while we mostly fix the value A e = 1 .1, we also vary the value between 0.8 and 1 .6 to see its systematic effect 
in §6. Incidentally this is useful in understanding the difference of the predicted number of arcs between open and Lambda-dominated 
CDM models found by Bartelmann et al. (1998). 



3. GRAVITATIONAL LENSING BY TRIAXIAL DARK MATTER HALOS 

In this section, we present several expressions for triaxial dark matter halos which are useful in calculating gravitational lensing 
properties. Under the thin lens approximation, gravitational lensing properties are fully characterized by the matter density projected 
along the line-of-sight (e.g., Schneider, Ehlers, & Falco 1992). We have to calculate the mass density profile projected along the 
arbitrary line-of-sight directions, because the line-of-sight, in general, does not coincide with the principal axis of a triaxial dark 
matter halo. 

For simplicity, in this section we redefine x/Rq and x'/Ro as x and x', respectively. In this case, 

<Wcrit(z) 



with R being defined by equation (4). In terms of the observer's coordinates (x',y',z'), R is written as 



R=Vfz' 2 +gz' + h, 



where 



/ = sin 2 



-rCOS 



j+T^sin 
b 



2 6] +cos 2 6> 



g = sin 9 sin 2<fi 



b 2 



h = ( — sin cos 



x' + sin26> 1- 



■COS — 



b 2 



sin 2 c/> /, 



Defining two new variables z\ and ( 



+ sin20cos6* ( — — -^ )x'y' ' + 



g 2 



cos 



-rCOS 

a 1 



>+-prsrn 
b 



+ sin 



we rewrite equation (16) as 



R=\z' 2 + C 2 - 



Then the convergence k can be expressed as a function of (: 

R ro ° D ~ ro ° 



K - 



J-oo ^cnt J-oo V J 



^crit J-oo ^crit J-oo Vf 

3 This expression looks different from its counterpart (eq. [21]) of JS02 for two reasons 



,2 + C 2)^ W 



/gNFw(C); 



(15) 

(16) 

(17) 
(18) 
(19) 

(20) 
(21) 

(22) 
(23) 



One is due to a typo in JS02 who omitted the factor 



y/ A v i r (z c ;JS02)/A v i r (z;JS02). Since A v ; r (JS02) = f2(z)A v j r according to the notation of this paper, this recovers the difference in the latter part. The other is 
the fact that we also incorporate the additional axis ratio dependence of c e which is noted in equation (23) of JS02. This explains the prefactor before A e in equation 
(14) of this paper. 
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where 

, _ 1 45 ce p C n t (z)Ro 

OTNFW = -J= =; , (24) 

V J ^cnt 

and 

,oo 1 

/gnfwW = / — — — — dz. (25) 

Jo (^Te) (i + \/^+?J 



The critical surface mass density E cr j t is defined by 

c 2 D 05 



^-AnGDoM (26) 

where Dol, Ads, and Dls denote the angular diameter distances from the observer to the lens plane, from the observer to the source 
plane, and from the lens plane to the source plane, respectively. 

The meaning of the variable ( can be easily understood by substituting equations (17)-(19) into equation (21): 

t 2 = j(Ax' 2 + Bx>y' + Cy> 2 ), (27) 



where 



2 2 \ 2 2 

C . 9 C i \ c c 



A = cos 2 9 sin 2 (j)+-r cos 2 (j> | + - T -rsin 2 6>, (28) 
b z J a 2 b 



2 2 
C L c 



B = cos0sin20 -,--5- , (29) 



2 2 

-r sin 2 (b+^r cos 2 6. (30) 
b 2 a 



7 2 

The quadratic form of equation (27) implies that the iso-£ curves are ellipses, and that the position angle of ellipses tp is 

1 B 

^=-arctan . (31) 

V 2 A-C 

By rotating the x'y'-plane by the angle tp, we diagonalize equation (27) such that 

r' 2 v' 2 

t 2 = ~ + y - (32) 

<& 1y 

q\ = 2/ = = , (33) 

A + C-^(A-C) 2 +B 2 

2 2/ 

ql = ; (34) 

y A + C+^/(A-C) 2 +B 2 



where 



Note that q x > q y for the given tp. We further define the axis ratio q as 

_ q y _ (A + C-^(A-C) 2 + B 2 Y /2 



(35) 



qx \A + C+^f(A-C) 2 + B 2 ) 

which represents the ellipticities of the projected isodensity curves of the triaxial dark halos. In this case, the convergence k is 
expressed as k = where £ 2 =x' 2 +y' 2 /q 2 . The advantage of this diagonalization is that we can apply the previous method 
to calculate lensing properties (Schramm 1990; Keeton 2001a) where the deflection angle (3 = (/5y,/3/) is expressed as a one- 
dimensional integral of the convergence k(£): 

PAx',y')=qx%(x',y'), (36) 
(3y{x',y')=qy'h{x',y'), (37) 

where the integral J n (x,y) is 

, f 1 «(g(v)) 

Jo [\-(l- q 2 )v ] 

and £(v) is 

e(v) = v(x 2 + / V (39) 

Figure 2 plots PDFs of q, q x , and q y . They were computed numerically using equations (1 1)-(12) and (33)-(35) under the assumption 
that the triaxial halo orientations (i.e., the angles 6 and (j>) are randomly distributed. In this plot we set M v i r = lO 15 /z -1 M and z = 0.3, 
which are s typical mass scale and a redshift of lensing clusters. It is clear from Figure 2 that the axis ratio of projected isodensity 
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0.2 0.4 0.6 0.8 1 

q 

FIG. 2. — PDFs of q (eq. [35]) , q x (eq. [33]), and q y (eq. [34]). Here we consider a halo with mass M v j r = 10 15 /r' Mq and redshift z = 0.3, but the result only 
weakly depends on the halo mass and redshift. These PDFs are calculated from PDFs of axis ratios p(a/c) and p(a/b) for which we use equations (11) and (12). 
We assume that the orientations of dark halos are random. 



contours strongly deviates from unity, having maximum around q ~ 0.6. This large degree of ellipticity suggests that the triaxial dark 
halos in realistic cosmological models significantly enhances the number of arcs compared with the conventional spherical model 
predictions. 

For a = 1, /gnfw(/) defined in equation (25) is analytically expressed as (Bartelmann 1996): 



/gnfwM : 



1 



l-r 2 
1 



■1 



-1 + - 



VT- 



zarctanh 



l-r 
1+7 



1 



: arc tan 



r-1 
7+1 



(r<D, 



(r>l), 



(40) 



but it does not has a simple analytical expression for a = 1 .5. Thus we use the following fitting formula in this case: 

2.614 



/gnfwM - 



r - 5 (l+2.378r 05833 + 2.617r 1 - 5 ) 



(41) 



The error of the above fit is < 0.6%. 



4. ARC STATISTICS IN THE TRIAXIAL DARK MATTER HALO 
4.1. Cross Sections for Arcs from the Monte Carlo Simulation 

First we compute the cross section for arcs without distinguishing tangential and radial arcs mainly because of the computational 
cost. In fact, the previous analyses indicate that while the number ratio of radial to tangential arcs offers another information on the 
density profile, the ratio is rather insensitive to the non-sphericity (Molikawa & Hattori 2001; Oguri et al. 2001; Oguri 2002a). 

Since the analytical computation of the cross sections is not practically feasible except for spherical models, we resort to the direct 
Monte Carlo method (Bartelmann & Weiss 1994; Miralda-Escude 1993b; Molikawa & Hattori 2001; Oguri 2002a). We showed that 
the convergence of triaxial dark matter halos is expressed by equation (23). Thus the corresponding lensing deflection angle a, and 
therefore the cross section a, are fully characterized by the two parameters, £>tnfw and q, as long as the finite size of source galaxies 
is safely neglected. Thus we perform the Monte Carlo simulations on the dimensionless X-Y plane, where X and Y are X = x' /(Roq x ) 
and Y = y' /(Roq x ), and tabulate the deflection angle and the dimensionless cross section 

a = a(b TNFW ,q), (42) 

a = CT(frrNFW,<7), (43) 

in 50 x 19 bins (a = 1) or 70 x 19 bins (a = 1.5) for £>tnfw and q, respectively. The dimensionless cross section is translated to the 
dimensional one in the source plane as 



a = a(b T NFw,q) x R q x 



DosV 
X D 0L 



(44) 



We follow the simulation method by Oguri (2002a) which is briefly summarized below. We use a 2048 x 2048 regular grid on 
the X-Y plane and calculate the deflection angle at each grid point. The box size is adjusted so as to include all arcs in the box for 
each (Z?tnfw,<7)- Therefore, the box size almost scales as the tangential critical line for each (£>tnfw 7 <?)- After those deflection angles 
are obtained at each grid point, we trace back the corresponding position in the source plane, and see whether or not it constitutes 
a part of lensed images. In order to take account of the source ellipticity which is also important in arc statistics (Keeton 2001b), 
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FIG. 3. — Contours of dimensionless cross sections <r and average magnification factors (fi) in the g-frpNFW plane for a = 1. The threshold axis ratios for arcs 
are set to e t i, = 10 (upper) and 1 {lower), respectively. Contours are drawn at 10 05 " for <r and at 10° 125,1 for (p), where n is integer. When n is in multiples of 4, 
contours are drawn by thick lines. These cross sections and magnification factors are derived from Monte Carlo simulations described in §4.1. 



we assume that it distributes randomly in the range of [0,0.5], where source ellipticity is defined by 1 -b s /a s with a s and b s being 
semi-major and semi-minor axes, respectively. We adopt this distribution of intrinsic ellipticities in order to compare our results with 
the previous works (e.g., Bartelmann et al. 1998) in which the same distribution was assumed. Moreover, the distribution is roughly 
consistent with the observed distribution (e.g., Lambas, Maddox, & Loveday 1992). Once we identify a lensed image, we compute 
its length I and width w as described in Oguri (2002a). Finally we define a lensed arc if the ratio of / and w exceeds the threshold 
value e,h that we set: 

- > eth- (45) 
w 

In practice, we consider e± = 7 and 10 to check the robustness of the conclusion. We also compute the average magnification of 
the arcs (/i) for each set of (£>tnfw,?) which is required in estimating the magnification bias (Turner 1980; Turner, Ostriker, & Gott 
1984). The contours of the lensing cross sections and the average magnification are plotted in Figures 3 and 4 for a = 1 and a = 1.5, 
respectively. We confirmed that the cross sections for q = 1 cases reproduce the analytic result of spherical lens models for point 
source. 

We should note that our current method does not take account of the finite size effect of source galaxies, and thus our results 
are, strictly speaking, applicable only to a sufficiently small source. Since the number of tangential arcs, which dominates the total 
number of arcs, is known to be insensitive to the source size (Hattori et al. 1997b; Bartelmann et al. 1998; Molikawa & Hattori 2001; 
Oguri et al. 2001; Oguri 2002a), this should not change our conclusion. 

4.2. Predicting Numbers of Arcs 

The next step is to average the cross section for arcs corresponding to a halo of M v ; r at zl and a galaxy at zs over the halo properties 
(its orientations and axis ratios): 



o-(M vil -,z L ,zs) = J d(a/c) J dc e J d(a/b) J dQ J d(j) p(a/c)p(c e )p(a/b\a/c)p(9)p((f>)a. 

In what follows, we assume the orientations of triaxial dark matter halos are completely random: 

sin 9 



P(0)- 
p{<t>)- 



2 

1 

2tt' 



(46) 

(47) 
(48) 
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FIG. 4. — Same as Figure 3, except for a = 1.5. 
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The realistic prediction for the number of arcs also requires to properly take account of the magnification bias. Thus we use the 
average of the cross section times number density of galaxies above the magnitude limit: 

OT%(Afvk,ZL,Zs) = J d(a/c) J dc e J d(ajb) J dO J d(j>p{a/c)p(c e )p{a/b\a/c)p(9)p((j))a J dLn g (L,zs), (49) 

where n„(L,z) is the luminosity function of source galaxies for which we adopt the Schechter form: 



n„(L,z)dL = 4>* 



exp 



L\dL 
Z* / Z*~' 



Its integral over L simply reduces to 



/•OO 

/ dLn g (L,zs) = 0*r(a s + l,L min /L*), 



(50) 



(51) 



with T(a,x) being the incomplete gamma function of the second kind. The lower limit of the integral, L\\ m , may be computed from 
limiting magnitude of observation, to*, and the lensing magnification factor (//) (see Figs. 3 and 4): 



^min 

~L*~ 



Y Q - 0.4(mii m — m* ) 



to* =M*+51og 



£>os(l+zs) 
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+ K(z s ). 



We adopt the K-correction in B-band for spiral galaxies (King & Ellis 1985): 



#(z) = -0.05 + 2.35z + 2.55z 2 -4.89z 3 +1.85z 4 . 

Finally the number distribution of lensed arcs for a halo of mass M v ; r at zl is given by 

dNmc cdt , 

— — (zs;M vir ,z L ) = crn g (M v ir,ZL,zs)^— (1+zs) , 
dzs ' dzs 



and the total number of lensed arcs for the halo is 

Narc(M vir ,Z L ): 



dzs —, — (zs;M vir ,z L )- 
dzs 



(52) 
(53) 
(54) 
(55) 

(56) 



While the upper limit of redshifts of source galaxies, zs,max, is in principle arbitrary, it is practically limited by the validity of the 
input luminosity function of source galaxies and the applied K-correction at high redshifts. In the present analysis, we conservatively 
set zs.max = 1-25 because of the K-correction (eq. [54]) and the luminosity function (§4.3). Nevertheless we stress here that our 
methodology can be applied to at higher redshifts if they are replaced by any reliable models valid there. 




4.3. Luminosity Function of Source Galaxies 

While the predicted number of arcs sensitively depends on the luminosity function of source galaxies (e.g., Hamana & Futamase 
1997), n g {L.z) is still fairly uncertain especially at high z. Thus we consider the following four luminosity functions measured up 
to z = 1.25: HDF1 from the Hubble Deep Field and the New Technology Telescope Deep Field (Poli et al. 2001), HDF2 from the 
Hubble Deep Field (Sawicki, Lin, & Yee 1997), SDF from the Subaru Deep Field (Kashikawa et al. 2003), and CFRS from the 
Canada-France Redshift Survey (Lilly et al. 1995). They are summarized in Table 1. Although the Schechter fits to those luminosity 
functions are valid only at z > (0.2 ~ 0.6), we simply extrapolate the values even down to z = if necessary. This does not affect our 
result in §5 at all since galaxies at z ~ 1 are the main sources of lensed arcs for our sample of clusters at z > 0.2 (§5.1). 

Except for HDF1, the Schechter parameters were derived assuming the Einstein-de Sitter (EdS) model (f2o = 1, Ao = 0) in the 
original references. We convert them into the counterparts in the Lambda-dominated universe (f^o = 0.3, Ao = 0.7) as follows. 

Since the number of galaxies in the redshift interval [zs,zs + dzs], 

2 cdt 

dN g (zs) oc D os -—dzsn g (L,zs)dL, (57) 
dzs 

is observable, it should be invariant. Thus the luminosity function in the Lambda-dominated universe is related to that in the EdS as: 

[n g (L',z s )dL'l = [D 2 0S (cdt/dz s )\ Ed5 r (L )dL , 5g) 

L g JLamMa ^OsM'AM Lambda ^ 



where 



[Dl 



[^OSjEdS 

The resulting luminosity functions in terms of the absolute magnitude M: 



jj _ l OS J Lambda (59) 



(M,z)rfM = O.9210*lO-°- 4(Qs+1)(M - M * ) exp -lO- - 4 ^*) ^, (60) 



at z = 0.3, 0.6, 0.9, and 1.2 are plotted in Figure 5. Clearly the uncertainty increases at fainter luminosities at z > 1, which may 
significantly change the predicted number of arcs. Therefore, while we adopt HDF1 as our fiducial model, we also attempt to 
evaluate the uncertainty due to the different choice of luminosity functions using the other three. 

4.4. Predicted Cross Sections and Numbers of Arcs 

Figure 6 shows the average cross sections (eq. [46]) of a dark matter halo of M vlI = 10 15 !i~ 1 Mq at zl = 0.3. The cross section for 
the triaxial model is larger by a factor of 10 (a = 1) and of 4 (a = 1 .5) than that for the spherical counterpart. Since the magnification 





factor is always larger for smaller cross sections (see Figures 3 and 4), the magnification bias further reduces the difference between 
a = 1 and 1 .5 for the triaxial model. This explains the behavior of Figure 7 where the source redshift distribution of arcs (eq. [55]) is 
plotted. Actually the figure indicates that the non-spherical effect even exceeds that of the difference due to the inner slope. 

Figures 8, 9 and 10 show how the predicted number of arcs depends on the mass of a lensing halo, the limiting magnitude of the 
survey, and the adopted luminosity function of source galaxies. Figure 8 shows that the number of arcs is sensitive to the mass of 
halo, implying the estimate of the mass of the target cluster is essential in interpreting the data. In addition, the difference between 
a = 1 and 1.5 becomes smaller for the triaxial model of M v ; r > 10 15 M Q . Thus in order to distinguish the inner slope clearly as well, 
one needs a sample of less massive clusters that have lensed arcs. 

Figure 9 indicates that the number of arcs is also sensitive to the magnitude limit, suggesting that the well-controlled selection 
function for the arc survey is quite important. On the other hand, the uncertainty of the luminosity function of source galaxies seems 
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FIG. 8. — Predicted numbers of arcs (eq. [56]) as a function of halo mass M yK . The redshift of the dark halo is still fixed to zl = 0.3. The B-band magnitude limit 
for arcs is set to m\\ m = 24. 
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FIG. 9. — Predicted numbers of arcs as a function of B-band magnitude limit m\\ m . The mass of lens cluster is M ylI = 10 l5 /i *Mq . 



to be less critical, at least for arcs of galaxies at z s < 1.25 that we consider in this paper (Fig. 10). The difference among the four 
luminosity functions (see Table 1) is merely up to 50 % for m\ lm < 24, and is within a factor of 2 even at my lm < 26 except CFRS. The 
predictions based on HDF1 approximately correspond to the median among the four and this is why we choose this as our fiducial 
model in what follows. 

5. COMPARISON WITH THE OBSERVED NUMBER OF ARCS 
5.1. Cluster Data 

We use a sample of 38 X-ray selected clusters compiled by Luppino et al. (1999). The clusters are selected from the Einstein 
Observatory Extended Medium Sensitivity Survey (EMSS). For all the clusters, deep imaging observations with B-band limiting 
magnitude m\ lm ~ 26.0 were carried out to search for arcs. 

As we remarked in the previous section, the mass estimate of those clusters is important in understanding the implications from 
the observed arcs statistics. For this purpose, we first construct a gas temperature - X-ray luminosity (in the Einstein band) relation 
from a subset of the above clusters whose temperature is determined. Then we estimate the temperature of the remaining clusters 
using the temperature - luminosity relation. Finally we estimate the mass of each cluster employing the virial mass - gas temperature 
relation of Finoguenov, Reiprich, & Bohringer (2001). 

More specifically, our best-fit luminosity - temperature relation from Figure 1 1 is 

'L z (0.3-3.5keV)\ 7 
lO^ergs" 1 ) ' 



ix.o 



(61) 



where 7 = 0.381 ±0.052 and Tx.o = 3.52^|gkeV. The derived luminosity - temperature relation is consistent with recent other 
estimations (e.g., Ikebe et al. 2002). Neglecting the possible redshift evolution for the luminosity - temperature relation (e.g., 
Mushotzky & Scharf 1997), we estimate the temperature of those clusters without spectroscopic data as shown in Table 2. The mass 
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FIG. 10. — Predicted numbers of arcs for different luminosity functions of source galaxies as a function of B-band magnitude limit mn m . Parameters of luminosity 
functions are given in Table 1, Only triaxial dark matter halo model and threshold axis ratio e t h = 10 are considered. 
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FIG. 1 1. — The luminosity - temperature relation for the EMSS cluster sample. Among 38 clusters, we use 21 clusters with measured temperature to derive 
luminosity - temperature relation. The best- fit luminosity - temperature relation is shown in equation (61). 



- temperature relation that we adopt is 

/ M \ 0,54 

This relation is derived by Shimizu et al. (2003) who converted the result of Finoguenov et al 
density profile (eq. [3]; the difference between a = 1 and 1.5 turned out to be negligible). 

5.2. Observed Number of A res 

The observed giant arcs (e t h = 10) in the 38 EMSS cluster sample are listed in Table 3. The number of arcs in this sample is roughly 
consistent with more recent data from different cluster samples (Zaritsky & Gonzalez 2003; Gladders et al. 2003). In order to be 
consistent with our adopted luminosity functions and K-correction of source galaxies, we need to select the arcs with z < 1.25. In 
reality, this is quite difficult; most of the observed arcs do not have a measured redshift, while uncertainties of source redshifts may 
systematically change lensing probabilities. For instance, Wambsganss et al. (2003) explicitly showed that it is important to take 
correctly account of the source redshift which can change cross sections by an order of magnitude. Moreover four in the list labeled 
"Candidate" in Table 3 are even controversial and may not be real lensed arcs. Thus we consider the two extreme cases; one is to 
select only the two arcs with measured redshifts less than 1.25, and the other is to assume that all the arcs without measured redshifts 
in the list (including the candidates) are located at z < 1.25. Of course the reality should be somewhere in between, and thus we 
assume that the range between the two cases well represents the current observational error. This means that the observational error 
can be greatly reduced if redshifts of all arcs are measured in the future observations. 



(62) 

. (2001) in terms of M v [ r assuming the 
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5.3. Comparison of Theoretical P redictions with Observations 

Finally let us compare our theoretical predictions with the data in detail. Our prediction of the number of arcs is the sum of 
equation (56) over all the 38 EMSS clusters: 

38 

Nemss = ^AWMot^zl,;)- (63) 
i=i 

We also compute the error of the predicted number of arcs by propagating the mass uncertainty for each cluster (Table 2). Figure 
12 shows the number of arcs in the 38 EMSS cluster sample as a function of the B-band limiting magnitude m\ lm . When the B-band 
magnitude of an arc is not available, we convert its corresponding V- or R-band magnitude into the B-band assuming typical colors 
of spiral galaxies at z ~ 1, B—V = V-R = 1 (Fukugita et al. 1995). 

The important conclusion that we draw from Figure 12 is that the triaxial model in the Lambda-dominated CDM universe with 
the inner slope of a = 1.5 successfully reproduces the observed number of arcs, and that the spherical model prediction with a = 1 
fails by a wide margin. Both the triaxial model with a = 1 and the spherical model with a = 1 .5 are marginal in a sense that the 
presence of substructure in the dark halo which we ignore in the current method should systematically increase our predicted number 
of arcs. Indeed Meneghetti, Bartelmann, & Moscardini (2003a) reported that the substructure enhances the number of arcs with 
e t h = 10 typically by a factor 2 or 3. This is exactly the amount of enhancement that is required to reconcile those two models with 
the observation. 

We note here that the additional contribution due to galaxies inside a cluster is generally small; Flores, Mailer, & Primack (2000) 
and Meneghetti et al. (2000) found that galaxies increase the number of arcs merely by ^10%. Even a central cD galaxy produces 
the number of arcs by not more than ~50% (Meneghetti, Bartelmann, & Moscardini 2003b). 

6. DISCUSSION 

6.1. Comparison with the previous result 

Our result that the halos in a Lambda-dominated CDM universe reproduces the observed number of arcs seems inconsistent with 
the previous result of Bartelmann et al. (1998) who claimed that only open CDM models can reproduce the observation. One 
possibility to explain the apparent discrepancy is the difference of the inner profile of halos; we showed that the slope of a = 1.5 
is required to reproduce the observation. This implies that N-body simulations may underestimate the real number of arcs unless 
they have sufficient spatial resolution. On the other hand, cluster-scale halos may indeed have a shallower inner profile (Jing & Suto 
2000). Therefore this is closely related to the well known problem of the inner slope of CDM dark matter halos (Navarro et al. 
1996, 1997; Moore et al. 1999; Jing & Suto 2000; Fukushige & Makino 2001, 2003; Power et al. 2003), and would need further 
investigation. Moreover, in reality, the mass estimate for each cluster, the limiting magnitude of source galaxies, and the adopted 
luminosity function would also affect the prediction in a more complicated fashion, and the further quantitative comparison is not 
easy at this point. 

Nevertheless we can point out the general tendency that open CDM models produce more arcs than Lambda-dominated CDM 
models because of the larger value of the concentration parameter in the former. Thus it is unlikely that difference between open and 
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FIG. 13. — The number of arcs in the 38 EMSS cluster sample as a function of A e for mi; m = 24. Dash-dotted line indicates the fiducial value for A e , A e = 1 . 1 , 
in a Lambda-dominated CDM model. Dotted lines suggest possible range of A e with taking account of the enhancement of the concentration parameter in an open 
CDM model (see text for details). 
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FIG. 14. — The number of arcs in the 38 EMSS cluster sample for fixed axis ratios of dark matter halos. The B-band limiting magnitude is set to m\ im = 24. Left 
panel plots the oblate case (a < b = c) while right panel is the prolate case (a = b < c). 



Lambda-dominated CDM models results from the "global" effect of the cosmological parameters. In order to show this, we compute 
the number of arcs as a function of A e still assuming the Lambda-dominated CDM model. Figure 13 plots A^emss for mn m = 24 as 
a function of A e . While JS02 found A e = 1.1 in a Lambda-dominated CDM models, their fitting formula (see also Bartelmann et 
al. 1998) tend to predict ~ 30-40% larger concentration parameter in open CDM models. This enhancement of the concentration 
parameter corresponds to A e = 1.43 ~ 1.53 if we still assume Lambda-dominated CDM models as a background cosmology. Thus 
the effect of A e alone increases the number of arc by ~ 50%- 100% even for triaxial cases, which is qualitatively consistent with the 
result of Bartelmann et al. (1998). 

6.2. Required Non-sphericity ofLensing Halos 

Although we showed that the triaxial halos predicted in the Lambda-dominated CDM model reproduce the observed number of 
arcs, the analysis employed a series of fairly complicated PDFs for the axial ratios of JS02, and it is not so clear what degree of 
non-sphericity for lensing halos is required to account for the observation. Thus we rather simplify the situation and consider that 
all halos consist of oblate (a < b = c) or prolate (a = b < c) halos with a fixed axial ratio. This is equivalent to replacing p(a/c) (eq. 
[11]) or p(a/b) (eq. [12]) by the corresponding ^-functions. Figure 14 plots the result of this exercise. 

The predicted number of arcs is indeed sensitive to the axis ratios of dark matter halos, and prolate halos of a/c < 0.5 in the a = 1 .5 
case reproduce the observation. This is basically consistent with the finding of JS02 for halo properties. 




10" 1 1 10 1 
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FIG. 15. — The effect of the sample variance. The number of arcs in the 38 EMSS cluster sample is calculated, not being averaged over axis ratios, orientation, 
and concentration parameters, but using fixed values for each clusters. Axis ratios, orientation, and concentration parameters for each cluster are randomly chosen 
according to their corresponding PDFs (§2.2). We calculate 10000 realizations and plot the histogram of frequency. Averaged values are shown by arrows. 



The reason why prolate halos tend to produce the larger number of arcs than oblate halos is explained as follows. Notice first 
that to keep the mass of dark matter halo invariant with the change of the axial ratio, £>tnfw should be approximately proportional 
to (ab/c 2 )~ l . Suppose that oblate and prolate halos are projected onto their axisymmetric direction (x for oblate and z for prolate). 
Then their lensing cross sections should scale as 

o-(oblate) oc CT((a/c)(a/c) _1 £>TNFW, 1) = o^tnfw, 1), (64) 

-J a((a/c)~ &tnfw,1)= y-J ctOtnfw,1), (65) 

where we assume <7(£>tnfW)<7) ^ ^tnfw- Since Figures 3 and 4 suggest 8 > 2, we find that cr(prolate) ^> er(oblate) for a/c < 1. If 
those halos are projected along the y-direction, on the other hand, their cross sections are almost the same: 

fa \ ~^ 

cr(oblate) ~ er(prolate) oc ^— J <7(Z?tnfw,<V c )- (66) 

The above consideration explains the qualitative difference between oblate and prolate halos, and points out that the elongation along 
the line-of-sight is also important in the arc statistics as well as the asymmetry of the projected mass density. 

6.3. Are Clusters Equilibrium Dark Matter Halos? 

So far we have assumed the one-to-one correspondence between dark matter halos and X-ray clusters. This assumption, however, 
is definitely over-simplified (Suto 2001, 2003). If "dark clusters" which are often reported from recent weak lensing analyses (Hattori 
et al. 1997a; Wittman et al. 2001; Miyazaki et al. 2002) are real, the one-to-one correspondence approximation may be unexpectedly 
inaccurate. As an extreme possibility, let us suppose that observed X-ray clusters preferentially correspond to halos in equilibrium. 
According to Jing (2000), such halos have generally larger concentration parameters and their scatter is small. In order to imitate 
this situation, we repeat the computation using A e =1.3 and the scatter of 0.18 (Jing 2000; JS02). We find that this modified model 
increases the number of arcs merely by 10%-20%. Thus our conclusion remains the same. 

6.4. Sample Variance 

The predicted number of arcs for the EMSS cluster that we have presented so far is based on the averaged cross section. This is 
a reasonably good approximation in the situation that the number of sample clusters is large enough, but in the current sample, its 
validity is not clear. To examine the sample variance, we re-compute the number of arcs in the 38 EMSS cluster sample without 
using the average statistics. Instead, we first randomly choose values of the axis ratios, the orientation angles, and the concentration 
parameters for each cluster according to their corresponding PDFs (§2.2). Then we sum up the number of arcs for the entire cluster 
sample. We repeat the procedure 10000 times each for a = 1 and 1.5, and construct a distribution function of A^emss as plotted in 
Figure 15. The resulting la sample variance is ^30% for a = 1 and ~15% for a = 1.5. Therefore we confirm that the effect of the 
sample variance does not change our overall conclusion. 

7. CONCLUSIONS 

We have presented a semi-analytic method to predict the number of lensed arcs, for the first time taking proper account of the 
triaxiality of lensing halos. We found that Lambda-dominated CDM models successfully reproduce the observed number of arcs 
of X-ray-selected clusters (Luppino et al. 1999) if the inner slope of the density profile is close to a = 1.5. Since the spherical 
models significantly underestimate the expected number of arcs, we conclude that the observed number of arcs indeed requires the 
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non-sphericity of the lensing halos. In fact, the number of arcs is sensitive to the axis ratios of those halos, and the non-sphericity that 
reproduces the observed number corresponds to the minor to major axis ratios of <~ 0.5. This value is perfectly consistent with the 
findings of JS02 in the Lambda-dominated CDM models. In this sense, we may even argue that the arc statistics lend strong support 
for the collisionless CDM paradigm at the mass scale of clusters. As discussed in Meneghetti et al. (2001), self-interacting dark 
matter models (Spergel & Steinhardt 2000) for instance, are inconsistent with the observed number of arcs not only because they 
erase the central cusp but because they produce much rounder dark matter halos (Yoshida et al. 2000a,b). Since we have exhibited 
that even the current arc surveys have a great impact in testing the collisionless CDM paradigm, larger surveys with well-controlled 
systematics in near future will unveil the nature of dark matter more precisely. 

Note added - In the published version, we used an incorrect PDF of the angle 9 (eq. [47]). In this version of the preprint, we used 
the correct PDF and replaced all related figures to the correct ones. In these new plots, however, the number of arcs in the triaxial 
dark matter halo model becomes smaller only by ~ 30%, so this does not affect the conclusion of the paper. 

We thank Y. P. Jing for useful correspondences concerning many aspects of the triaxial dark matter halo model from N-body sim- 
ulations, and Masahiro Takada for discussions. We also thank an anonymous referee for many useful comments. J. L. acknowledges 
gratefully the research grant of the JSPS (Japan Society of Promotion of Science) fellowship. This research was supported in part by 
the Grant-in-Aid for Scientific Research of JSPS (12640231). 
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ARC STATISTICS IN TRIAXIAL DARK MATTER HALOS 



Table 1 

b-band luminosity functions of source galaxies used in this paper. 



Name 


Model 


z Range 




M^- 5 log ft" 


0*[/z 3 Mpc- 3 ] 


Ref. 


HDF1 


Lambda b 


0.00 - 0.50 c 


-1.19 


-20.26 


2.5 x 10" 2 


1 






0.50 - 0.75 


-1.19 


-19.97 


2.9 x 10" 2 








0.75 - 1.25 


-1.25 


-20.61 


1.2 x 10" 2 




HDF2 


EdS d 


0.00 - 0.50 c 


-1.40 


-21.20 


9.0 x 10" 3 


2 






0.50- 1.00 


-1.30 


-19.90 


4.2 x 10" 2 








1.00- 1.25 


-1.60 


-22.10 


6.0 x 10" 3 




SDF 


EdS d 


0.00- 1.00 c 


-1.07 


-19.78 


4.2 x 10" 2 


3 






1.00- 1.25 


-0.92 


-20.13 


4.3 x 10" 2 




CFRS 


EdS d 


0.00 - 0.50 c 


-1.03 


-19.53 


2.7 x 10" 2 


4 






0.50 - 0.75 


-0.50 


-19.32 


6.2 x 10" 2 








0.75 - 1.00 


-1.28 


-19.73 


5.4 x 10" 2 








1.00- 1.25 f 


-2.50 


-21.36 


9.6 x 10" 4 





a B-band AB magnitude can be converted to conventional Johnson-Morgan mag- 
nitude via B AB = B-0.14 (Fukugita, Shimasaku, & Ichikawa 1995). 

b n Q = 0.3, A = 0.7. 
c Extrapolated to z = 0. 

d O = 1, A =0. 

e The luminosity function for blue galaxies only. 

References. — (1) Poli et al. 2001; (3) Sawicki et al. 1997; (3) Kashikawa et al. 
2003; (4) Lilly et al. 1995 
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Table 2 



Properties of clusters in the 38 EMSS distant cluster sample. 







L x (EdS) a 




Lx(Lambda) b 


T x 


Myir 




Name 




[10 44 ergs _1 ] 


[10 44 ergs _1 ] 


[keV] 


[10 14 /z _1 M Q ] 


Ref. 


MS 001 1 7+0837 

1V_I_VJ \J\J 11./ 1 V_/ 1_> / 


0.163 


3.77 


2.24 


4 79 + 2'l? c 


3 89 + 2 7 l 


1 


MS 001 5 9+1609 

iV_LO \J\J J. ^3 . 7T J. \J\jy 


0.546 


14.64 


11.03 


8 92 + S I 7 


12 30 + MS 


1, 2 


MS 0302 5+1717 


0.425 


2.88 


2.04 


4 62 + iM 5c 


3 64+2 

J -°^-0.60 




MS 0302 7+1658 


0.426 


5.04 


3.57 


4 Q<+0.80 
^■ JJ -0.64 


3 25+Hl 

-'• z ' J -0.89 


1, 2 


MS 0353.6-3642 


0.320 


5.24 


3.48 


6.46 + 'i^ 

-0.80 


6.77 + }?° 

— 1 55 


1, 2 


MS 0433 9+0957 


0.159 


4.34 


2.57 


5 O4. +0 ' 53c 
J - UH -0.48 


4 O7+0:83 




MS 0440 5+0204 


0.190 


4.01 


2.43 


5 30+2 -52 


4 69 +0 ' 98 


1, 3 


MS 0451 5+0250 


0.202 


6.98 


4.27 


8 60 + S S 

°- uu -0.50 


11 50 + Hl 


1, 3 


MS 0451.6-0305 


0.539 


19.98 


15.00 


-0.80 


15.97i :4 5 

-2.30 


1, 2 


MS 0735.6+7421 


0.216 


6.12 


3.79 


c oc+0.68c 

j,OJ r\ £-i 

-0.61 


"^-1.09 


1 


MS 081 1 6+6301 


0.312 


2.10 


1.40 


4 87 + 2'^ 

^■°'-0.63 


4 01 + i 


1, 2 


MS 0839 8+2938 


0.194 


5.35 


3.25 


^•^"-0.20 


3 05 + iH 7 


1, 3 


MS 0906 5+1 1 10 


0.180 


5.77 


3.47 


5 65 + 2 K c 


5 28 +L U 


1 


MS 1006.0+1201 


0.221 


4.82 


2.99 


5.34W 

-0.52 


-0.86 


1 


MS 1008 1-1224 


0.301 


4.49 


2.95 


8 21 +115 


10 55+H 4 . 


1, 2 


MS 1054 5-0321 


0.823 


9.28 


7.79 


10 4+1.00 
lu -^-i.oo 


16 35 +2 91 

iu.jj_2 


1, 4 


MS 1137.5+6625 


0.782 


7.56 


6.26 


5.70^ 

—0.60 


5.37 + [- 4 ° 

-1.05 


1, 5 


MS 1 147 3+1 103 


0.303 


2.30 


1.51 


5 96 +0 " 


5 83 +179 


1, 2 


MS 1201.5+2824 


0.167 


2.03 


1.21 


o to+0.34c 
-0.32 


9 s i +0.42 
-0.39 


1 


MS 1 208 7+3928 


0.340 


2.03 


1.37 


o Q7+0.36c 
J -^'-0.33 


o 7^+0.46 
z - ' J -0.42 




MS 1224 7+2007 


0.327 


4.61 


3.08 


4 09+2'S 
^• u ^-0.52 


2 90+2 


1, 2 


MS 1231.3+1542 


0.238 


2.88 


1.81 


4.41 + °- 42 ) c 

-0.39 


3.34+^ 9 

-0.55 


1 


MS 1241 5+1710 


0.549 


10.70 


8.07 


6 09 +L38 

°- u5, -1.14 


6 07 +2 ' 55 

°- u '-2.10 


1 2 


MS 1 244 2+71 14 


0.225 


3.84 


2.39 


4 QA+0.50c 


4 06 +077 

+ - uo -0.71 




MS 1253 9+0456 


0.230 


3.14 


1.96 


^■- ,J -0.40 


3 54 + S S 

J -'^-0.58 


1 


MS 1 358 4+6245 


0.327 


10.62 


7.09 


7 50 +4 ,'2 0d 


o Qo+9.48 
°-^ J -2.01 


1 2 6 


MS 1426 4+01 58 


0.320 


3.71 


2.47 


6 38+ ? 8 , 


6 62+H 8 . 

u ' uz, -2.30 


1, 2 


MS 1455.0+2232 


0.259 


16.03 


10.23 


— 1.15 


5.20 + ^3 

-1 .98 


1, 3,6 


MS 1 512 4+3647 


0.372 


4.81 


3.30 


3 39+R- 40 

- , --- ,5 '-0.35 


2 05 + 2'^ 

z - UJ -0.39 


1, 2 


MS 1546.8+1132 


0.226 


2.94 


1.83 


4 4 o+0.43c 
^■^ J -0.39 


a 07+0.61 
J - :, '-0.55 


1 


MS 1618.9+2552 


0.161 


2.24 


1.33 


o Q9+0.36c 
-'•^^-0.33 


2.68!^ 


1 


MS 1621.5+2640 


0.426 


4.55 


3.22 


6 59+ - 92 

°- ,y -0.81 


7.02!! : 8 g 


1,2 


MS 1910.5+6736 


0.246 


4.39 


2.78 




4 <Q+0.89 
+ - ,J -0.81 


1 


MS 2053.7-0449 


0.583 


5.78 


4.43 






1,2 


MS 2137.3-2353 


0.313 


15.62 


10.34 


j; on+1.09d 
J-^-0.42 


4 <a + 1.76 
^■ JJ -0.68 


1,2,6 


MS 2255.7+2039 


0.288 


2.04 


1.33 


o Q9+0.36c 
- , - : ^-0.33 


2-68!°;^ 


1 


MS 2301.3+1506 


0.247 


3.29 


2.08 


4 fi <r+0.46c 




1 


MS 2318.7-2328 


0.187 


6.84 


4.14 


fi c+0.73c 
°- UJ -0.65 


6.00!H 4 


1 



Note. — Errors are at 68% confidence limit. 

a X-ray luminosity in the 0.3 - 3.5keV band for Slo = 1, Ao = 0, and h = 0.5 universe. 
b X-ray luminosity in the 0.3 - 3.5keV band for fio = 0.3, Ao = 0.7, and h = 0.7 universe. 
c Estimated from Lx~Tx relation (eq. [61]). 
d The effects of cooling flows are corrected. 

References. — (1) Luppino et al. 1999; (2) Novicki, Sornig, & Henry 2002; (3) Mushotzky & 
Scharf 1997; (4) Jeltema et al. 2001; (5) Borgani et al. 2001; (6) Allen & Fabian 1998 



ARC STATISTICS IN TRIAXIAL DARK MATTER HALOS 



Table 3 

Giant arcs (l/w > 10) in the 38 EMSS distant cluster sample. 



Cluster 


Zh 


Arc 


zs 


l/w 




m 


arc 


Notes 


Ref. 


MS 0302.7+1658 


0.426 


Al 


~0.8 a 


> 18 


B 




23.8 




1,2 






A1W 




> 12 


B 




24.9 






MS 0440.5+0204 


0.190 


Al 


0.532 


> 10 


B 




22.9 




1, 3,4 






A3 




>20 


B 




24.0 






MS 0451.6-0305 


0.539 


Al 




10 


V 




24.6 




1 


MS 1006.0+1201 


0.221 


A2+A3 




>20 


V 


< 


22.1 


Candidate 


1,5 






A4 




12.9 


V 




21.4 


Candidate 




MS 1008.1-1224 


0.301 


A2 




10.0 


B 




23.4 


Candidate 


1 


MS 1358.4+6245 


0.328 


Al 


4.92 


>21 










1,6 


MS 1621.5+2640 


0.426 


Al 




> 18 


B 




23.1 




1,7 


MS 1910.5+6736 


0.246 


Al 




10.5 


R 




20.6 


Candidate 


1,5 


MS 2053.7-0449 


0.583 


AB 




>22 


V 




22.4 




1,7 


MS 2137.3-2353 


0.313 


Al 


1.501 


18.1 


B 




22.0 




1,8,9, 10 



'Estimated from color of the arc. 



References. — (1) Luppino et al. 1999; (2) Mathez et al. 1992; (3) Luppino et al. 1993; (4) 
Gioia et al. 1998 (5); Le Fevre et al. 1994; (6) Franx et al. 1997; (7) Luppino & Gioia 1992; (8) 
Fort et al. 1992; (9) Hammer et al. 1997; (10) Sand et al. 2002 



